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Abstract 

Gravitational radiation reaction forces and balance equations for energy and 
momenta are investigated to 3/2 post-Newtonian (1.5PN) order beyond the 
quadrupole approximation, corresponding to the 4PN order in the equations 
of motion of an isolated system. By matching a post-Newtonian solution for 
the gravitational field inside the system to a post-Minkowskian solution (ob- 
tained in a previous work (T|) for the gravitational field exterior to the system, 
we determine the 1PN relativistic corrections to the "Newtonian" radiation 
reaction potential of Burke and Thorne. The 1PN reaction potential involves 
both scalar and vectorial components, with the scalar component depending 
on the mass-type quadrupole and octupole moments of the system, and the 
vectorial component depending in particular on the current-type quadrupole 
moment. In the case of binary systems, the 1PN radiation reaction potential 
has been shown [QQ to yield consistent results for the 3.5PN approximation 
in the binary's equations of motion. Adding up the effects of tails, the radia- 
tion reaction is then written to 1.5PN order. In this paper, we establish the 
validity to 1.5PN order, for general systems, of the balance equations relating 
the losses of energy, linear momentum, and angular momentum in the system 
to the corresponding fluxes in the radiation field far from the system. 



I. INTRODUCTION 

The old idea (in any field theory) of losses of energy and momenta in an isolated system, 
due to the presence of radiation reaction forces in the equations of motion, is of topical 
interest in the case of the gravitational field. Notably, gravitational radiation reaction forces 
play an important role in astrophysical binary systems of compact objects (neutron stars or 
black holes). The electromagnetic-based observations of the Hulse-Taylor [|J binary pulsar 
PSR 1913+16 have yielded evidence that the binding energy of the pulsar and its companion 
decreases because of gravitational radiation reaction |j]-|8[]. 

Even more relevant to the problem of radiation reaction are the future gravitational- 
based observations of inspiralling (and then coalescing) compact binaries. The dynamics 
of these systems is entirely driven by gravitational radiation reaction forces. The future 
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detectors such as LIGO and VIRGO should observe the gravitational waves emitted during 
the terminal phase, starting about twenty thousands orbital rotations before the coalescence 
of two neutron stars. Because inspiralling compact binaries are very relativistic, and thanks 
to the large number of observed rotations, the output of the detectors should be compared 
with a very precise expectation from general relativity |9|-|TT||. In particular Cutler et al || 
have shown that our a priori knowledge of the relativistic (or post-Newtonian) corrections 
in the radiation reaction forces will play a crucial role in our ability to satisfactorily extract 
information from the gravitational signals. Basically, the reaction forces inflect the time 
evolution of the binary's orbital phase, which can be determined very precisely because of 
the accumulation of observed rotations. The theoretical problem of the phase evolution has 
been addressed using black-hole perturbation techniques, valid when the mass of one body 
is small as compared with the other mass [|12]-rr^|, and using the post-Newtonian theory, 



valid for arbitrary mass ratios fll8|- pi)1 . It has been shown []T3| , [T4| , |TT[| that post-Newtonian 
corrections in the radiation reaction forces should be known up to at least the third post- 
Newtonian (3PN) order, or relative order c~ 6 in the velocity of light. 

The radiation reaction forces in the equations of motion of a self-gravitating system arise 
at the 2.5PN order (or c -5 order) beyond the Newtonian acceleration. Controlling the nth 
post-Newtonian corrections in the reaction force means, therefore, controlling the (n+2.5)th 
post-Newtonian corrections in the equations of motion. If n = 3 this is very demanding, 
and beyond our present knowledge. A way out of this problem is to assume the validity 
of a balance equation for energy, which permits relating the mechanical energy loss in the 
system to the corresponding flux of radiation far from the system. Using such a balance 
equation necessitates the knowledge of the equations of motion up to the rzPN order instead 
of the (n + 2.5)PN one. The price to be paid for this saving is the computation of the far- 
zone flux up to the same (relative) nPN order. However, this is in general less demanding 
than going to (n + 2.5)PN order in the equations of motion. All the theoretical works on 
inspiralling binaries compute the phase evolution from the energy balance equation. Black- 
hole perturbations [|14}-|T6|| reach n = 4 in this way, and the post-Newtonian theory [|I~8| -|2T| 
has n = 2.5. 

An important theoretical problem is therefore to improve the present situation by showing 
the validity to post-Newtonian order of the balance equations for energy, and also for linear 
and angular momenta. This problem is equivalent to controlling the radiation reaction forces 
at the same post-Newtonian order. Arguably, this problem is also important in its own (not 
only for applications to inspiralling compact binaries). 

Radiation reaction forces in general relativity have long been investigated (see [22 for a 
review of works prior the seventies). In the late sixties, Burke and Thorne [P3]-[21)[|, using a 
method of matched asymptotic expansions, introduced a quasi-Newtonian reactive potential, 
proportional to the fifth time-derivative of the Newtonian quadrupole moment of the source. 
At about the same time, Chandrasekhar and collaborators [p^-|28|, pursuing a systematic 
post-Newtonian expansion in the case of extended fluid systems, found some reactive terms 
in the equations of motion at the 2.5PN approximation. The reactive forces are different in 
the two approaches because of the use of different coordinate systems, but both yield secular 
losses of energy and angular momentum in agreement with the standard Einstein quadrupole 
formulas (see however [[29|-|31"]1). These results, after later confirmation and improvements 
^|-PU[, show the validity of the balance equations to Newtonian order, and in the case 
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of weakly self-gravitating fluid systems. In the case of binary systems of compact objects 
(such as the binary pulsar and inspiralling binaries), the Newtonian balance equations are 
also known to be valid, as the complete dynamics of these systems has been worked out 
(by Damour and Deruelle |40|-(42|j22ll) up to the 2.5PN order where radiation reaction effects 
appear. 

Post-Newtonian corrections in the radiation reaction force can be obtained from first 
principles using a combination of analytic approximation methods. The methods are (i) a 
post-Minkowskian or non-linear expansion method for the field in the weak-field domain of 
the source (including the regions far from the source), (ii) a multipolar expansion method 
for each coefficient of the post-Minkowskian expansion in the domain exterior to the source, 
and (iii) a post-Newtonian expansion method (or expansion when c — ■> oo) in the near-zone 
of the source (including its interior). Then an asymptotic matching (in the spirit of |23|-p5|) 
permits us to connect the external field to the field inside the source. Notably, the methods 
(i) and (ii) have been developed by Blanchet and Damour p3|-fi5|| on foundations laid by 
Bonnor and collaborators f46-48||, and Thorne | |49|| . The method (iii) and matching have 
also been developed within the present approach |50|-p2 



The post-Newtonian correction that is due to gravitational wave tails in the reaction 
force was determined first using the latter methods [Q. The tails of waves are produced 
by scattering of the linear waves off the static spacetime curvature generated by the total 
mass of the source (see e.g. |47|j48[| ). Tails appear as non-local integrals, depending on the 



full past history of the system, and modifying its present dynamics by a post-Newtonian 
correction of 1.5PN order in the radiation reaction force, corresponding to 4PN order in the 
equations of motion [p3 |. It has been shown in that the tail contribution in the reaction 
force is such that the balance equation for energy is verified for this particular effect. This 
is a strong indication that the balance equations are actually valid beyond the Newtonian 
order (1.5PN order in this case). For completeness we shall include this result in the present 
paper. 

The methods (i)-(ii) have been implemented in [1] in order to investigate systematically 
the occurence and structure of the contributions in the exterior field which are expected 
to yield radiation reaction effects (after application of the method (iii) and the relevant 
matching). The present paper is the direct continuation of the paper [1], that we shall refer 
here to as paper I. 

Working first within the linearized theory, we investigated in paper I the "antisymmetric" 
component of the exterior field, a solution of the d'Alembertian equation composed of a 
retarded (multipolar) wave minus the corresponding advanced wave. Antisymmetric waves 
in the exterior field are expected to yield radiation reaction effects in the dynamics of the 
source. Indeed, these waves change sign when we reverse the condition of retarded potentials 
into the advanced condition (in the linearized theory), and have the property of being 
regular all over the source (when the radial coordinate r — > 0). Thus, by matching, the 
antisymmetric waves in the exterior field are necessarily present in the interior field as 
well, and can be interpreted as radiation reaction potentials. In a particular coordinate 
system suited to the (exterior) near zone of the source (and constructed in paper I), the 
antisymmetric waves define a radiation reaction tensor potential in the linearized theory, 



generalizing the radiation reaction scalar potential of Burke and Thorne [23 25 



Working to non-linear orders in the post-Minkowskian approximation, we introduced in 
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paper I a particular decomposition of the retarded integral into the sum of an "instanta- 
neous" integral, and an homogeneous solution composed of antisymmetric waves (in the 
same sense as in the linearized theory). The latter waves are associated with radiation 
reaction effects of non-linear origin. For instance, they contain the non-linear tail contri- 
bution obtained previously At the 1PN order, the non-linear effects lead simply to a 
re-definition of the multipole moments which parametrize the linearized radiation reaction 
potential. However, the radiation reaction potential at 1PN order has been derived only in 
the external field. Thus, it was emphasized in paper I that in order to meaningfully interpret 
the physical effects of radiation reaction, it is necessary to complete this derivation by an 
explicit matching to the field inside the source. 

We perform the relevant matching (at 1PN order) in the present paper. Namely, we 
obtain a solution for the field inside the source (satisfying the non- vacuum field equations), 
which can be transformed by means of a suitable coordinate transformation (in the exterior 
near-zone of the source) into the exterior field determined in paper I. The matching yields in 
particular the multipole moments parametrizing the reaction potential as explicit integrals 
over the matter fields in the source. As the exterior field satisfies physically sensible boundary 
conditions at infinity (viz the no-incoming radiation condition imposed at past-null infinity), 
the lPN-accurate radiation reaction potentials are, indeed, appropriate for the description 
of the dynamics of an isolated system. 

To 1PN order beyond the Burke-Thorne term, the reaction potential involves a scalar 
potential, depending on the mass-type quadrupole and octupole moments of the source, and 
a vectorial potential, depending in particular on the current-type quadrupole moment. The 
existence of such vectorial component was first noticed in the physically restricting case 
where the dominant quadrupolar radiation is suppressed [|T] . 



A different approach to the problem of radiation reaction has been proposed by Iyer and 
Will |fj,§] in the case of binary systems of point particles. The expression of the radiation 
reaction force is deduced, in this approach, from the assumption that the balance equations 
for energy and angular momentum are correct (the angular momentum balance equation 
being necessary for non-circular orbits). Iyer and Will determine in this way the 2.5PN and 
3.5PN approximations in the equations of motion of the binary, up to exactly the freedom 
left by the un-specified coordinate system. They also check that the lPN-accurate radiation 
reaction potentials of the present paper (and paper I) correspond in their formalism, when 
specialized to binary systems, to a unique and consistent choice of a coordinate system. This 
represents a non-trivial check of the validity of the 1PN reaction potentials. 

In the present paper we prove that the lPN-accurate radiation reaction force in the 
equations of motion of a general system extracts energy, linear momentum and angular 
momentum from the system at the same rate as given by the (known) formulas for the 
corresponding radiation fluxes at infinity. The result is extended to include the tails at 
1.5PN order. Thus we prove the validity up to the 1.5PN order of the energy and momenta 
balance equations (which were previously known to hold at Newtonian order, and for the 
specific effects of tails at 1.5PN order). 

Of particular interest is the loss of linear momentum, which can be viewed as a "recoil" 
of the center of mass of the source in reaction to the wave emission. This effect is purely 
due to the 1PN corrections in the radiation reaction potential, and notably to its vectorial 
component (the Newtonian reaction potential predicts no recoil). Numerous authors have 
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obtained this effect by computing the flux of linear momentum at infinity, and then by relying 
on the balance equation to get the actual recoil []54j -57]]. Peres |58|] made a direct computation 
of the linear momentum loss in the source, but limited to the case of the linearized theory. 
Here we prove the balance equation for linear momentum in the full non-linear theory. 

The results of this paper apply to a weakly self-gravitating system. The case of a source 
made of strongly self-gravitating (compact) objects is a priori excluded. However, the 
theoretical works on the Newtonian radiation reaction in the binary system PSR 1913+16 
40|- [42tf22f have shown that some "effacement" of the internal structure of the compact bodies 



is at work in general relativity. Furthermore, the computation of the radiation reaction at 
1PN order in the case of two point-masses has shown agreement with a formal reduction, 
by means of ^-functions, of the 1PN radiation reaction potentials, initially derived in this 
paper only in the case of weakly self-gravitating systems. These works give us hope that 
the results of this paper will remain unchanged in the case of systems containing compact 
objects. If this is the case, the present derivation of the 1.5PN balance equations constitutes 
a clear support of the usual way of computing the orbital phase evolution of inspiralling 
compact binaries [p|-[2T|. 

The plan of this paper is the following. The next section (II) is devoted to several recalls 
from paper I which are necessary in order that the present paper be essentially self-contained. 
In Section III we obtain, using the matching procedure, the gravitational field inside the 
source, including the 1PN reactive contributions. Finally, in Section IV, we show that the 
latter reactive contributions, when substituted into the local equations of motion of the 
source, yield the expected 1PN and then 1.5PN balance equations for energy and momenta. 



II. TIME-ASYMMETRIC STRUCTURE OF GRAVITATIONAL RADIATION 

A. Antisymmetric waves in the linearized metric 

Let D e = {(x, t), |x| > r e } be the domain exterior to the source, defined by r e > a, where 
a is the radius of the source. We assume that the gravitational field is weak everywhere, 
inside and outside the source. In particular a ^> GM/c 2 , where M is the total mass of 
the source. Let us consider, in D e , the gravitational field at the first approximation in a 
non-linearity expansion. We write the components h^ v of the deviation of the metric density 
from the Minkowski metric 7]^ u in the form 

bT = V^gg^ - jf v = Gh% + 0{G 2 ) , (2.1) 



'(i) 

where the coefficient of the Newton constant G represents the linearized field , satisfying 
the vacuum linearized field equations in D e , 

Dhft = d»d x h% + d v d x h x ^ - TTdxdah^ . (2.2) 

We denote by □ = rj^ u d^d u the flat space-time d'Alembertian operator. 

The general solution of the equations ( |2.2|) in D e can be parametrized (modulo an ar- 
bitrary linearized coordinate transformation) by means of two and only two sets of multi- 
pole moments, referred to as the mass- type moments, denoted Ml, and the current-type 



moments, Sl H49fl . The capital letter L represents a multi-index composed of I indices, 
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L = %\%% ■ ■ -i\ (see |59| for our notation and conventions). The multipolarity of the moments 
is / > in the case of the mass moments, and Z > 1 in the case of the current moments. 
The M^'s and S^s are symmetric and trace- free (STF) with respect to their I indices. The 
lowest-order moments M, Mi and Si are constant, and equal respectively to the total con- 
stant mass (including the energy of the radiation to be emitted), to the position of the center 
of mass times the mass, and to the total angular momentum of the source. The higher-order 
moments, having I > 2, are arbitrary functions of time, Mi(t) and Sl(£), which encode all 
the physical properties of the source as seen in the exterior (linearized) field. In terms of 
these multipole moments, the "canonical" linearized solution of Thorne [Ml reads 



ft can(l)- --^l^—n-OL 
u l>0 
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(2.3c) 



This solution satisfies ( |2.2|) and the condition of harmonic coordinates (i.e. d u h^ n ^ = 0). 
Here we impose that the multipole moments M^(t) and «Sz,(t) are constant in the remote past, 
before some finite instant — T in the past. With this assumption the linearized field fl2.3| ) 
(and all the subsequent non-linear iterations built on it) is stationary in a neigthborhood 
of past-null infinity and spatial infinity (it satisfies time-asymmetric boundary conditions in 
space-time). This ensures that there is no radiation incoming on the system, which would 
be produced by some sources located at infinity. 

In paper I (Ref. |l|]) the contribution in ( |2.3|) which changes sign when we reverse the 
condition of retarded potentials to the advanced condition was investigated. This contribu- 
tion is obtained by replacing each retarded wave in ( |2.3| ) by the corresponding antisymmetric 
wave, half difference between the retarded wave and the corresponding advanced one. The 
antisymmetric wave changes sign when we reverse the time evolution of the moments Mi(t) 
and Siit), say Mi(t) — > Ml(— t), and evaluate afterwards the wave at the reversed time —t. 
Thus, (|2.3j) is decomposed as 



The symmetric part is given by 

(^can(l)) 



antisym 



( h 0i 



' sym ' — ' Z' 

l>0 



4 ^ tl dL f M dt ~ r /c) + M L(t + r/c) 
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I 4 v (-)'Z ( ^ L _ 1 (t-r/c) + ^ L . 1 (t + r/c) l 
+ ^^(7TT)!^ aL - 1 l 2r / ' (2 ' 5b) 

f hij \ _^_(-V \ Mlf L _ 2 {t-r/c) + Mif L _ 2 {t + r/c) 

-^ 2 w^ daL - 2 x aKi 2^ / • (2 - 5c) 

The antisymmetric part is given similarly. However, as show in paper I, it can be re-written 
profitably in the equivalent form 

The second, third and fourth terms clearly represent a linear gauge transformation, asso- 
ciated with the gauge vector . This vector is made of antisymmetric waves, and reads 

* -c^^TKT3T)^{ £ j ' (2 ' 7a) 

^(-)' (2Z + l)(2i + 3) a fi\4- 2) (t-r/c)-Mi- 2) (t + r/c) 
^ = ~ 2 "li 77/ — T\ diL 



l>2 



l\ 1(1 - 1) 2r 



, 4 (-)' 2/ + 1 / M iL „x(t - r/c) - M <L _i(t + r/c) 1 

+ ? g "TTT^T j 27 | 

+ rr 8 ^- 1 j 27 j • (2 - 7b) 

Note that even though we have introduced first and second time-antiderivatives of the mul- 
tipole moments, denoted e.g. by m£ (t) and m| ^(t), the dependence of ( |2.7|) on the 
multipole moments ranges in fact only in the time interval between t — r/c and t + r/c 
(see paper I). The first term in ( |2.6| ) defines, for our purpose, a radiation reaction tensor 
potential VJ^ C in the linearized theory (in this term s takes the values 0, 1,2 according to 
\iv = 00, Oi, ij). The components of this potential are given by [^J 

reac " h~ ] - — i('-i) I 27 / ' (2 - 8a) 

T/ o* ^yH' (/ + 2)(2/ + l) a (^(t-r/cJ-^-^t+r/c) 
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(see Eqs. (2.19) of paper I). By adding the contributions of the gauge terms associated with 
( |2.7|) to the radiation reaction potential ( |2.8|) one reconstructs precisely, as stated by Q2.6|) , 
the antisymmetric part (/i^nOantisym °f the linearized field ( |2.3| ). 

The scalar, vector, and tensor components of (|2.8| ) generalize, within the linearized the- 
ory, the Burke-Thorne p3|-p5|1 scalar potential by taking into account all multipolarities of 
waves, and, in principle, all orders in the post-Newtonian expansion. Actually a full justifica- 
tion of this assertion would necessitate a matching to the field inside the source, such as the 
one we perform in this paper at 1PN order. At the "Newtonian" order, the 00 component 
of the potential reduces to the Burke-Thorne potential, 



00 
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(t)+0 



(2.9) 



At this order the Oi and ij components make negligible contributions. Recall that a well- 
known property of the Burke-Thorne reactive potential is to yield an energy loss in agreement 
with the Einstein quadrupole formula, even though it is derived, in this particular coordinate 



system, within the linearized theory (see [30|). In this paper we shall show that the same 
property remains essentially true at the 1PN order. [This property is in general false for other 
reactive potentials, valid in other coordinate systems, for which the non-linear contributions 
play an important role.] Evaluating the reaction potential VJ^. at the 1PN order beyond 
051), we find that both the 00 and Oi components are to be considered, and are given by 
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45 



(2.10a) 
(2.10b) 



At this order the ij components of the potential can be neglected. 

In the next subsection we address the question of the corrections to the 1PN reaction 
potential fl2.10p which arise from the non-linear contributions to the exterior field. Answering 
this question means controlling the non-linear metric at the 3.5PN order. 



B. The 3.5PN approximation in the exterior metric 

The radiation reaction potential ( p.8| )-( |2TT0| ) represents the antisymmetric part of the 
linearized metric in a particular coordinate system, obtained from the initial harmonic co- 
ordinate system in which Q2.3p holds by applying the gauge transformation associated with 
(|2.7|) . In this coordinate system, the new linearized metric reads 

Ki) = KU) + <9T + d»e - rr^e . (2.H) 

It fulfills, of course, the linearized equations (|2.2|) . Furthermore, since = 0, the harmonic 
coordinate condition is still satisfied. However, we recall from paper I that the latter new 
harmonic coordinate system is not completely satisfying for general purposes, and should 
be replaced by a certain modified (non-harmonic) coordinate system. The reason is that the 
gauge vector defined by ( |2.7|) is made of antisymmetric waves, and consequently the metric 
( |2-11| ) contains both retarded and advanced waves. In particular, the metric is no longer 
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stationary in the remote past (before the instant — T where the moments are assumed to 
be constant). Of course, the advanced waves which have been introduced are pure gauge. 
Nevertheless, the non-stationarity of the linearized metric in the remote past breaks one 
of our initial assumptions, and this can be a source of problems when performing the non- 
linear iterations of the metric by this method. Therefore, it was found necessary in paper I to 
replace the gauge vector fl2.7Q by a modified gauge vector, such that the modified coordinate 
system has two properties. First, it reduces in the near- zone of the source (the domain 
Di defined in Sect. Ill), to the un-modified coordinate system given by (|2.7| ), with a given 
but arbitrary post-Newtonian precision. Second, it reduces exactly to the initial harmonic 
coordinate system in which the linearized metric is ( p.3[ ) in a domain exterior to some time- 
like world tube surrounding the source (in fact, a future-oriented time-like cone whose vertex 
is at the event t = —T, x = 0). This modified coordinate system is non-harmonic. It has 
been suggested to the author by T. Damour in a private communication, and is defined in 
Sect. II. C of paper I. By the first property, we see that if we are interested to the field in 
the near- zone of the source, and work at some finite post-Newtonian order (like the 1PN 
order investigated in this paper), we can make all computations using the un-modified gauge 
vector (p.7|) . Indeed, it suffices to adjust a certain constant, denoted K in paper I, so that the 
modified gauge vector agrees with ( |2.7| ) with a higher post-Newtonian precision. Thus, in 
the present paper, we shall not use explicitly the modified coordinate system. By the second 
property, we see that the standard fall-off behavior of the metric at the various infinities 
(notably the standard no-incoming radiation condition at past-null infinity) are preserved 
in the modified coordinate system. By these two properties, one can argue that what only 
matters is the existence of such a modified coordinate system, which permits us to make the 
connection with the field at infinity, but that in all practical computations of the metric in 
the near-zone one can use the un-modified coordinate system defined by (2.7). 

Based on the linearized metric fl2.11| ) (or, rather, on the modified linearized metric (2.29) 
of paper I) we built a full non-linear expansion, 

hT = GTig + G 2 h^ + G 3 h^ + ... , (2.12) 

satisfying the vacuum field equations in a perturbative sense (equating term by term the 
coefficients of equal powers of G in both sides of the equations). The non-linear coefficients 
^(2)> ^(3)' ••• are n ^ e ^(T) * n ^ ne form of multipole expansions parametrized by and 
Sl- The construction of the non-linear metric is based on the method of ||43|| . As we are 



considering multipole expansions valid in D e (and singular at the spatial origin r = 0), we 
need to use at each non-linear iterations of the field equations a special operator generalizing 
the usual retarded integral operator when acting on multipole expansions. We denote this 
operator by JFD^ 1 , to mean the "Finite part of the retarded integral operator" (see |43|] for 
its precise definition). The non- linear coefficients hf%>, WZ, ... are given by 

h% = F n^A^h^) + q$ , (2.13a) 
h^=^a R 'A^(h {1) ,h i2) ) + q^ , ... , (2.13b) 

where the non-linear source terms A^j, Aj^, ... represent the field non-linearities in vacuum, 
and depend, at each non-linear order, on the coefficients of the previous orders. The second 
terms q^, q^, ■■■ ensure the satisfaction of the harmonic coordinate condition at each 
non- linear order (see [[43|1 ). 



9 



When investigating the 3.5PN approximation, we can disregard purely non-linear effects, 
such as the tail effect, which give irreducibly non-local contributions in the metric inside the 
source. These effects arise at the 4PN approximation (see |)3 and Sect. IV below). Still 
there are some non- linear contributions in the metric at the 3.5PN approximation, which 
are contained in the first two non-linear coefficients h%. and h/%\ given by (|2.13| ). These 
contributions involve some non-local integrals, but which ultimately do not enter the inner 
metric (after matching). As shown in paper I, the contributions due to and in the 
1PN radiation reaction potential imply only a modification of the multipole moments Ml 
and Sl parametrizing the potential. We define two new sets of multipole moments, 



M L {t) = M L {t) + 



S L (t) = S L (t) + 



£?m(t) for / = 
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where the functions m, mi and s< are given by the non-local expressions 
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(3), 



+lf dv f dw\-^M^(w)M^\w)-% iab M^(w)S l 

& J-oo J-oo L 63 45 



w: 



i(t) = -\e iab f dv M£{v)M h { f{v) + Hi(t) . 

D J — oo 



(2.15a) 



c z 



(2.15b) 
(2.15c) 



The function T L {t) in (|2.14a|) , and the functions F(t), Gi(t) and Hi(t) in ( p.!5|) , are some 
local (or instantaneous) functions, which do not play a very important role physically (they 
are computed in paper I). Then the radiation reaction potential at the 1PN order, in the 
non- linear theory, is given by the same expression as in ( [2.10P , but expressed in terms of the 
new multipole moments, say M. = {Ml,Sl}, 



V- c [M] = -^x a x b MH\t) + ^ 



+0 



1 x abc M$ c (t)-±-r*x ab MZ\t) 



189 



vr c \M\ = % 



^ iah M§\t)-^e i ^Sg\t) 



70 



+ 04 



(2.16a) 
(2.16b) 



(see Eq. (3.53) in paper I). In the considered coordinate system, the metric, accurate to 
1PN order as concerns both the usual non-radiative effects and the radiation reaction effects, 
reads (coming back to the usual covariant metric g e ^) 
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rfff = -i + l (y cxt [M] + V vcac [M]) - 1 (y cxt [M] + V Icac [M] 

+^*& + ^*& + o(± 5 ) , (2.17a) 

4 (vrm + vr^iM}) + ± ^ + 4 rc 4 + (4) > ( 2 - 17b ) 
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o ext = 5- ■ 



1 + 4 (v ext [A4] + F reac [A4] 



The superscript ext is to remember that the metric is valid in the exterior domain D e , 
and will differ from the inner metric by a coordinate transformation (see Sect. III). The 
Newtonian and 1PN approximations are entirely contained in the externals potentials V cxt 
and V^ ext , given by the multipole expansions of symmetric waves, 



\M L (t- r -)+M L (t , , . 
V^[M] = GJ2 y -^-d L | —± cJ 2r °J- \ , (2.18a) 




_ r_y f m$ x ft - r ) + i ft + 1 



■^W^ 1 { | • (2 ' 18b) 

The 2PN and 3PN approximations are not controlled at this stage; they are symbolized 
in ( j2.17|) by the terms c~ n n g Qx ^. However, these approximations are non-radiative (or non- 
dissipative), as are the Newtonian and 1PN approximations (the 1PN, 2PN, and 3PN terms 
are "even" in the sense that they yield only even powers of 1/c in the equations of motion). 
A discussion of the 4PN approximation in the exterior metric can be found in Sect. III.D of 
paper I. 

As it stands, the metric ( |2.17| ) is disconnected from the actual source of radiation. The 



multipole moments Ml and Sl are left as some un-specified functions of time. Therefore, 
in order to determine the radiation reaction potentials ( |2.16| ) as some explicit functional of 
the matter variables, we need to relate the exterior metric ( |2.17 ) to a metric valid inside 



the source, solution of the non-vacuum Einstein field equations. We perform the relevant 
computation in the next section, and obtain the multipole moments Ml and Sl as integrals 
over the source. 



III. THE 1PN-ACCURATE RADIATION REACTION POTENTIALS 

A. The inner gravitational field 

The near-zone of the source is defined in the usual way as being an inner domain Di = 
{(x, t), |x| < n}, whose radius satisfies r» > a (Di covers entirely the source), and r« <C X 
(the domain Di is of small extent as compared with one wavelength of the radiation). These 
two demands are possible simultaneously when the source is slowly moving, i.e. when there 
exists a small parameter of the order of 1/c when c — > oo, in which case we can assume 



11 



Ti/X = 0(1/ c). Furthermore, we can adjust r e and r« so that a < r e < r, (where r e defines 
the external domain D e in Sect. II). 

In this subsection we present the result of the expression of the metric in Di, to the 
relevant approximation. In the next subsection we prove that this metric matches to the 
exterior metric reviewed in Sect. II. The accuracy of the inner metric is 1PN for the usual non- 
radiative approximations, and 1PN for the dominant radiation reaction. Thus, the metric 
enables one to control, in the equations of motion, the Newtonian acceleration followed by 
the first relativistic correction, which is of order c~ 2 or 1PN, then the dominant "Newtonian" 
radiation reaction, of order c~ 5 or 2.5PN, and finally the first relativistic 1PN correction in 
the reaction, c~ 7 or 3.5PN. The intermediate approximations c~ 4 and c -6 (2PN and 3PN) 
are left un-determined, like in the exterior metric (2.17). The inner metric, in Di, reads 



flSS = -1 + |v m - |(V in ) 2 + i flfl S + ± a «fc + O (i) , (3.1a) 
& = + + ^ + O (i) , (3.1b) 



It is valid in a particular Cartesian coordinate system (x, t), which is to be determined by 
matching. Like in (2.17), the terms c~ n n g l // u represent the 2PN and 3PN approximations. 
Note that these terms depend functionally on the source's variables through some spatial 
integrals, extending over the whole three-dimensional space, but that they do not involve any 



non-local integral in time. These terms are "instantaneous" (in the terminology of |53fl) and 
"even" , so they remain invariant in a time reversal, and do not yield any radiation reaction 
effects. The remainder terms in ( |3.1| ) represent the 4PN and higher approximations. Note 
also that some logarithms of c arise starting at the 4PN approximation. For simplicity we do 
not indicate in the remainders the dependence on In c. The potentials V m and Vf 1 introduced 
in ( |3.1| ) are given, like in fl2.17| ), as the linear combination of two types of potentials. With 
the notation V™ = (V m ,Vj- n ), where the index /x takes the values 0,i and where = V in 
we have 



V; = V^K] + V^ C [1] . (3.2) 

The first type of potential, V^ n , is given by an integral of the symmetric potentials, i.e. 
by the half-sum of the retarded integral and of the corresponding advanced integral. Our 
terminology, which means here something slightly different from Sect. II, should be clear 
from the context. We are referring here to the formal structure of the integral, made of the 
sum of the retarded and advanced integrals. However, the real behavior of the symmetric 
integral under a time-reversal operation may be more complicated than a simple invariance. 
The mass and current densities = (a, a/) of the source are defined by 

^OCI _|_ rpkk 

cr = 5 , (3.3a) 

cr 

o-j = — , (3.3b) 
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where T^ v denotes the usual stress-energy tensor of the matter fields (with T the spatial 
trace £<5j-fcT J ). The powers of 1/c in ( |3.3|) are such that a M admits a finite non-zero limit 



when c — > +00. The potentials V? n are given by 



d 3 x' 



<7„ 



x , i |x 



x', i H — |x 

c 



(3.4) 



To lowest order when c — > +00, y m reduces to the usual Newtonian potential, and VJ m to the 



usual gravitomagnetic potential. It was noticed in |50| that when using the mass density a 
given by (3.3a), the first (non-radiative) post-Newtonian approximation takes a very simple 
form, involving simply the square of the potential in the 00 component of the metric. See 
also (4.5) below, where we use the post-Newtonian expansion V m = U + d 2 X/2c 2 + 0(c -4 ). 

The fact that the inner metric contains some symmetric integrals, and therefore some 
advanced integrals, does not mean that the field violates the condition of retarded potentials. 
The metric (371) is in the form of a post-Newtonian expansion, which is valid only in the 
near zone Z)j. It is well-known that the coefficients of the powers of 1/c in a post-Newtonian 
expansion typically diverge at spatial infinity. This is no concern of us because the expansion 
is not valid at infinity (it would give poor results when compared to an exact solution). Thus, 
the symmetric integral ( p.4| ) should more properly be replaced by its formal post-Newtonian 
expansion, readily obtained by expanding by means of Taylor's formula the retarded and 
advanced arguments when c — ► +00. Denoting d 2p = (d/dt) 2p we have 



V m ( 



+00 1 

x, t) = G Y . . , „ / d 3 x'|x - x'P-^VJx', t) . 



(3.5) 



This expansion could be limited to the precision indicated in ( |3.1|) . It involves (explicitly) 
only even powers of c -1 , and is thus expected to yield essentially non-dissipative effects. 
However, the dependence of V^ n on c -1 is more complicated than indicated in (|3.5|) . Indeed, 
the mass and current densities depend on the metric (|3.1|), and thus depend on c _1 
starting at the post-Newtonian level. Even more, the densities and their time-derivatives 
do contain, through the contribution of the reactive potentials V^ eac (see below), some odd 
powers of c" 1 which are associated (a priori) to radiation reaction effects. These "odd" 
contributions in V^ n form an integral part of the equations of motion of binary systems at 
the 3.5PN approximation ||, but we shall prove in Sect. IV that they do not participate to 
the losses of energy and momenta by radiation at the 1PN order, as they enter the balance 
equations only in the form of some total time-derivatives. The secular losses of energy and 
momenta are driven by the radiation reaction potentials V7^ eac , to which we now turn. 

The lPN-accurate reaction potentials V^ eac = (y reac , V* ea,c ) involve dominantly some odd 
powers of c _1 , which correspond in the metric (|3.1|) to the 2.5PN and 3.5PN approximations 
c~ 5 and c~ 7 taking place between the (non-dissipative) 2PN and 3PN approximations. Since 
the reactive potentials are added linearly to the potentials V^ n , the simple form mentionned 
above of the 1PN non-radiative approximation holds also for the 1PN radiative approxi- 
mation, in this coordinate system. The y^ eac 's are given by exactly the same expressions 
as obtained in paper I for the exterior metric (see ( |2.16| ) in Sect. II), but they depend on 
some specific "source" multipole moments I = {II, Jl} instead of the unknown multipole 
moments A4. Namely, 
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Oli 



17' 



x >*) = -r 



01 %ijkljk(t) A r £ ijk x jmJkm(f) 



(5), 



21 



45 



(3.6a) 
(3.6b) 



where we recall our notation x^k = x^ — -^{SijXf. + Sik%j + ^jk x i) [ 59fl . The multipole 
moments Iij(t), Iijk(t) and Jjj(t) are some explicit functionals of the densities er M . Only 
the mass quadrupole ijj(t) in the first term of V reac needs to be given at 1PN order. The 
relevant expression is 



1 2~ o2 20 



d x { XijU H -x Xijd+ a — 

1 « 14c 2 »j * 



21c 2 



Xijkdt<Jk 



(3.7) 



(see ( |3.21a| ) for the general expression of II). The mass octupole and current quadrupole 
lij kit) and Jijif) take their standard Newtonian expressions, 



hjk = J d 3 yix ijk cr + O (—^j 
d 3 xe kr 



J, 



km<i x j>k@m 



(3.8a) 
(3.8b) 



The potentials V rea,c and V^ reac generalize to 1PN order the scalar reactive potential of Burke 
and Thorne [p^ -25], whose form is that of the first term in ( |3.6a ). The vectorial potential 
-y^reac g^gj-g ^ e equations of motion at the same 3.5PN order as the 1PN corrections in 
V reac . The first term which is neglected in V rea,c , of order c~ 8 or 1.5PN, is due to the tails 
of waves (see Sect. IV. C). 



B. Matching to the exterior field 

In this subsection we prove that the inner metric presented above (i) satisfies the Einstein 
field equations within the source (in the near- zone _Dj), and (ii) matches to the exterior metric 
( [2.17] ) in the intersecting region between and the exterior zone D e (exterior near-zone 

An£> e ). 

Note that during Proof (i) we do not check any boundary conditions satisfied by the 
metric at infinity. Simply we prove that the metric satisfies the field equations term by term 
in the post-Newtonian expansion, but at this stage the metric could be made of a mixture 
of retarded and advanced solutions. Only during Proof (ii) does one check that the metric 
comes from the re-expansion when c — > oo of a solution of the Einstein field equations sat- 
isfying some relevant time-asymmetric boundary conditions at infinity. Indeed, the exterior 
metric has been constructed in paper I by means of a post-Minkowskian algorithm valid 
all over the exterior region D e , and having a no-incoming radiation condition built into it 
(indeed, the exterior metric was assumed to be stationary in the remote past, see Sect. II). 

The proof that ( |3.1| ) is a solution admissible in follows immediately from the particular 
form taken by the Einstein field equations when developed to 1PN order [|50| 

□ ln(-^ n ) = + O evcn (1) , (3.9a) 
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□flS = + Oevcn (i) , (3.9b) 

□fl& = -^?V + Oe W mf^ • (3.9C) 



The point is that with the introduction of the logarithm of — as a new variable in (|3.9aj) 



the equations at the 1PN order take the form of linear wave equations. The other point is 
that the neglected post-Newtonian terms in ( |3.9|) are "even", in the sense that the explicit 
powers of c _1 they contain, which come from the differentiations of the metric with respect to 
the time coordinate x° = ct, correspond formally to integer post-Newtonian approximations 
(to remember this we have added the subscript "even" on the O-symbols). This feature is 
simply the consequence of the time-symmetry of the field equations, implying that to each 
solution of the equations is associated another solution obtained from it by a time-reversal. 

Because the potentials V™ satisfy exactly DV^ 11 = — AnGa^, a consistent solution of ( p.9| ) 
is easily seen to be given by Q3.1|) in which the reactive potentials V^ eac are set to zero. Now 
the equations ( |3.9j ) are linear wave equations, so we can add linearly to V^ n any homogeneous 
solution of the wave equation which is regular in Di. One can check from their definition 
( WM that the reactive potentials V^ eac form such a homogeneous solution, as they satisfy 

n yre,c = q f_\ ^ (3 1Qa ) 

□V7 eac = O (^j . (3.10b) 

So we can add V^ eac to V^ n , defining an equally consistent solution of (|3.9|), which is precisely 
Q3.1|) , modulo the error terms coming from ( [3.10D and which correspond to the neglected 



4PN approximation. [Note that V^ eac comes from the expansion of the tensor potential 
of the linearized theory, which satisfies exactly the source-free wave equation. It would be 
possible to define V^ eac in such a way that there are no error terms in ( |3.10| ). See (|4.33| ) for 
more precise expressions of V^ eac , satisfying more precisely the wave equation.] 

With Proof (i) done, we undertake Proof (ii). More precisely, we show that ( |3.1| ) differs 
from the exterior metric ( |2.17| ) by a mere coordinate transformation in D,ir\D e . This will be 
true only if the multipole moments Ml and Sl parametrizing ( |2.17| ) agree, to the relevant 



order, with some source multipole moments II and Jl- Fulfilling these matching conditions 
will ensure (in this approximate framework) the existence and consistency of a solution of 
the field equations valid everywhere in Di and D e . As recalled in the introduction, this 
is part of the method to work out first the exterior metric leaving the multipole moments 
arbitrary (paper I), and then to obtain by matching the expressions of these moments as 
integrals over the source (this paper). 

To implement the matching we expand the inner metric ( |3.1| ) into multipole moments 
outside the compact support of the source. The comparison can then be made with ( |2.17|) , 



which is already in the form of a multipole expansion. Only the potentials V^ n need to 
be expanded into multipoles, as the reactive potentials V^ eac are already in the required 
form. The multipole expansion of the retarded integral of a compact-supported source is 
well-known. E.g., the formula has been obtained in the Appendix B of |50| using the STF 



formalism for spherical harmonics. The multipole expansion corresponding to an advanced 
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integral follows simply from the replacement c — > — c in the formula. The script letter Ai 
will be used to denote the multipole expansion. M.{V™) reads as 



M(V [n )=Gj2 [ - 1 rd L 



l>0 



M(Vh=Gj2 y -^d L 



l>0 



' F L (t-r/c)+F L (t + r/c) \ 
2r J 

'G iL (t-r/c) + G iL (t + r/c) 
2r 



(3.11a) 
(3.11b) 



where Fi(t) and Gn{t) are some tensorial functions of time given by the integrals 

Flit) = J d 3 x.XL J dz 5i(z)a(x, t + z\x\/c) , 
Git(t) — J d 3 xx L J ^dzSi(z)ai(x,t + z\x\/c) . 



(3.12a) 
(3.12b) 



The function 8i(z) appearing here takes into account the delays in the propagation of the 
waves inside the source. It reads 



W = ^f(i-* 2 )'; fjz5 l (z) = i 



(3.13) 



(see Eq. (B.12) in |50|). Note that the same functions F^t) and Gn(t) parametrize both the 
retarded and the corresponding advanced waves in ( ^.11| ). Indeed 5i (z) is an even function 
of its variable z, so the integrals (|3.12|) are invariant under the replacement c — > —c. 

Using an approach similar to the one employed in |50 , we perform an irreducible decom- 
position of the tensorial function Gn (which is STF with respect to its I indices L but not 
with respect to its Z + 1 indices iL), as a sum of STF tensors of multipolarities I + 1, I and 
I — 1. The equation (2.17) in |50[ gives this decomposition as 



G 



iL 



I 21- 1 



l + l 



21 + 1 



(3.14) 



where the tensors Cl+i, Dl and El_\ (which are STF with respect to all their indices) are 
given by 



G L +i{t) = J d 3 x J dz6i(z)x <L (Ti l+1> {x,t + z\x\/c) , 

D L{t) = J d 3 x J ^dz5i(z)e ab<il x L ^ 1>a a b {x.,t + z\x\/c) 
E L -i(t) = J c? 3 x J ^dzdi(z)x aL -i(Ta(x,t + z\x\/c) . 



Then by introducing the new definitions of STF tensors 

A L = F L 
B L = IC L - 



4(2/ + 1) (i) 
c 2 (/ + l)(2/ + 3) L ' 



I 



c 2 (/ + l)(2/ + 3) 



E {2) 

L > 



(3.15a) 
(3.15b) 
(3.15c) 

(3.16a) 
(3.16b) 
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and by using standard manipulations on STF tensors, we can re-write the multipole expan- 
sions (|3.11| ) in the new form 

M(V-) = -eW + CE^j Mt = '' /C) + } ■ (3.17a) 

Miyh = ^ _ G gV s " { B ' x - i(< "' /c) 2 + r B,L " (<+r/c) } 

_ G g ii^,, { ^ - :M±£^ ± ^ ) , (3 . 17b ) 

where we denote 

„ 4G^ (-)' 21 + 1 \ E L (t-rlc) + Edt + rlc) \ 

<t> = ~ g (myraTs 91 j S ) ' (318) 

Next the moment is expanded when c — > oo to 1PN order, and the moments B^, to 
Newtonian order. The required formula is (B.14) in [ 50|1 , which immediately gives 

A fj |- 1 2~ o2 4(2/ +1) A 1 

Ai = 7 d X \ XLU + 2c°(2i + 3) X ^ ° ~ c^ + l)(2l + 3) X ^] 

+O even (^) , (3.19a) 

B L = l J rf 3 x£ <L _KT i;> + O even , (3.19b) 

D l = J d 3 yLE ab<il x L ^ 1>a (T b + O even (—^j . (3.19c) 

Here the notation O even (c _n ) for the post-Newtonian remainders simply indicates that the 
whole post-Newtonian expansion is composed only of even powers of c _1 , like in ( |3.5|) (the 
source densities <t m being considered to be independent of c" 1 ), as clear from (B.14) in pD| . 
We now transform the leading-order term in the equation for Bl using the equation of 
continuity for the mass density a. The Newtonian equation of continuity does the needed 
transformation, but one must be careful about the higher-order post-Newtonian corrections 
which involve some reactive contributions. It can be checked that these reactive contributions 
arise only at the order 0(c~ 7 ), so that the equation of continuity reads, with evident notation, 
d t a + di<Ji = O even (c~ 2 ) + 0(c~ 7 ). From this one deduces 



d 

Bl = dt 



| J d 3 xx L (x j + O even (i) + O Qf\ . (3.20) 

All elements are now in hands in order to compare, in the exterior near- zone Di D D e , 
the metrics ( |3.1| ) and ( |2.17| ). The "source" multipole moments X = {II, Jl} are defined by 
the dominant terms in ( |3.19a|) and ( |3.19c|) , 

f f 1 4(2/ + 1) ] 

I L = j rf 3 x |x L a + 2c2(2j + 3) x 2 ^V - c2{l + 1){2l + 3 f iLd t «j , (3.21a) 

J L = / d 3 xe ab<il x L ^ 1>a cr b . (3.21b) 
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The mass-type moment It includes 1PN corrections, while the current-type moment Jt is 



Newtonian. The mass moment It was obtained in ||50|| , where it was shown to parametrize 
the asymptotic metric generated by the source at the 1PN order. When I = 2 and I = 3 we 
recover the moments introduced in (|3.7|) - (|3.8|) . With (|3.19|) - (|3.21|) , the multipole expansions 
( p. 17] ) become 



l>0 



2r 



4^ "^"F^l 2r 

n ^{-Y l a f JbL-i(t-r/c) + J bL -i(t + r/c)\ 

— 2^ ,, , , ; £iabOaL-l 



2r 



+0, 



even i n 



(3.22a) 



(3.22b) 



Thus, from the definition ( 2.18Q of the external potentials V^ xt , we obtain the relationships 



M(V m ) = -cd t <f + V cxt [l] + o cvcn ( I 



from which we readily infer that the multipole expansion M.(g™ v ) of the metric 
in Di n D P , 



(3.23a) 
(3.23b) 
If) reads, 



M(g*) + ~d t < 



M(g$)+d i( p 



-l + ^(v ext [i] + v vcac \i]) - ^(v cxt [i] + v icac [i]y 



1 _ in 1 _ in 



l + -(^ ext [X] + \/ rcac [J]) 



+ + ^ + o n 



(3.24a) 
(3.24b) 
(3.24c) 



Clearly the terms depending on 0° have the form of an infinitesimal gauge transformation 
of the time coordinate. We can check that the corresponding coordinate transformation can 
be treated, to the considered order, in a linearized way (recall from ( |3.18|) that 4>° is of order 
c~ 3 ). Finally, in the "exterior" coordinates 



X 



cxt 



X i + O evcn -T J + O (s 



(3.25a) 
(3.25b) 



the metric ( |3.24| ) is transformed into the "exterior" metric 
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g*f = -l + -(v cxt [i] + v Tcac [i}) - -(v cxt [i] + v reac [x]Y 

+^ + ^ + °(^o) , (3.26a) 

<7 e f = -4(v; cxt [x] + vr c [i\) + 4r f + + o (i) , (3.26b) 



o cxt = 5 
iJii u 



1.1 



1 + ^(V cxt [l] + V Icac [l]) 



+ + ^ + 0(1) . (3.26c-) 



This metric is exactly identical, as concerns the 1PN, 2.5PN, and 3.5PN approximations, to 
the exterior metric ( [2.17|) obtained in paper I, except that here the metric is parametrized by 
the known multipole moments X instead of the arbitrary moments Ai. Thus, we conclude 
that the two metrics ( |3.1| ) and ( p. 17] ) match in the overlapping region Di H D e if (and only 



if) there is agreement between both types of multipole moments. This determines Ml and 
Sl- More precisely, we find that M L and Sl must be related to II and Jl given in ( |3.21| ) 
by 

M L = I L + O cvcn + O , (3.27a) 

S L = J L + O even (\)+o(\) , (3.27b) 



,c 2 J Vc 6 , 

where as usual the relation for Ml is accurate to 1PN order, and the relation for Sl is 
Newtonian (we do also control the parity of some neglected terms). Satisfying the latter 
matching solves the problem at hand, by showing that the inner metric (|3.1|) - (|3.8|) results 
from the post-Newtonian expansion of a solution of the (non-linear) field equations and a 
condition of no incoming radiation. 

We emphasize the dependence of the result on the coordinate system. Of course, the 
metric (|3.1| ), which contains the reactive potentials ( |3.6| )- (|3.8| ), is valid only in its own 
coordinate system. It is a well-known consequence of the equivalence principle that radiation 
reaction forces in general relativity are inherently dependent on the coordinate system (see 
e.g. for a comparison between various expressions of the radiation reaction force at the 
Newtonian order). The coordinate system in which the reactive potentials ( |3.6j )- (|3lf ) are 
valid is defined as follows. We start from the particular coordinate system in which the 
linearized metric reads ( |2.3j ). Then we apply two successive coordinate transformations. 
The first one is associated with the gauge vector £ M given by (|2.7f ), and the second one 
is associated with whose only needed component is 0° given by ( |3.18| ). The resulting 
coordinate system is the one in which V^ eac is valid. (Actually, the gauge vector £ M should 
be modified according to the procedure defined in Sect. II. C of paper I, so that the good 
fall-off properties of the metric at infinity are preserved.) 

IV. THE BALANCE EQUATIONS TO POST-NEWTONIAN ORDER 

A. Conservation laws for energy and momenta at 1PN order 

Up to the second post-Newtonian approximation of general relativity, an isolated system 
admits some conserved energy, linear momentum, and angular momentum. These have been 
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obtained, in the case of weakly self-gravitating fluid systems, by Chandrasekhar and Nutku 



27 1 . The less accurate lPN-conserved quantities were obtained before, notably by Fock |6T 



In this subsection we re-derive, within the present framework (using in particular the mass 
density a defined in ( |3.3a|) ), the lPN-conserved energy and momenta of the system. The 
1PN energy and momenta are needed in the next subsection, in which we establish their 
laws of variation during the emission of radiation at 1PN order (hence we do not need the 
more occurate 2PN-conserved quantities). 
To 1PN order the metric (13. II) reduces to 
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where V m and V™ are given by ( |3.4j ). In fact, V m and V™ are given by their post-Newtonian 
expansions ( |3.5|) , which can be limited here to the terms 
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where the instantaneous (Poisson-like) potentials U, X and U{ are defined by 
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Since V m is a symmetric integral, there are no terms of order c -3 in ( |4.2a| ) (such a term 
would be a simple function of time in the case of a retarded integral). We shall need (only 
in this subsection) a metric whose space-space components ij are more accurate than in 
( 4.1c ), taking into account the next-order correction term. We introduce an instantaneous 
potential whose source is the sum of the matter stresses, say = T % \ and the (Newtonian) 
gravitational stresses, 
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The metric which is accurate enough for our purpose reads, in terms of the instantaneous 
potentials flOft-flPD, 
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9oo = -1 + + ±[d?X - 2U 2 ] + O (A) , (4.5a) 
C = -|^ + 0(^) , (4.5b) 
flS = (l + f ^ + ^[^^ + 2 ^]) + - SijP] + O (1) . (4.5c) 



The square-root of (minus) the determinant of the metric is 

2- l^,,on2 ' 



.gin = 1 + _ U + _ [d f X + 2U 2 - AP] + O (^J . (4.5d) 

Consider the local equations of motion of the source, which state the conservation in the 
covariant sense of the stress-energy tensor T^ u (i.e. V M X^ = 0). These equations, written 
in a form adequate for our purpose, are 

d^ a = T a , (4.6) 

where the left-hand-side is the divergence in the ordinary sense of the material stress-energy 
density 



and where the right-hand-side can be viewed as the four-force density 

1 
2 



T a = -yj-^T^d a ^ (4.8) 



The lPN-conserved energy and momenta follow from integration of these equations over the 
ordinary three-dimensional space, which yields the following three laws (using the Gauss 
theorem to discard some divergences of compact-supported terms) 

j t [- J Anj} =-cj d 3 xf , (4.9a) 

ay***} (49b » 

(1 -e ijk [ d 3 x Xj U k ) =e ijk [ d 3 x (xjP k + U{) . (4.9c) 



dt [c 

The quantities IT^ and T a are then determined. With (|4.5|) we obtain, for the various 
components of IT^, 

U° = -ac 2 + an + ^[ a P-aiUi]+o(^\ , (4.10a) 

n° = c^(l+4^)-^ + o(i) , (4.10b) 

n o = (l " 4^) - + O (I) , (4.10c) 

Hj. = a l3 (l + 4 £) - ±a,f/, + O (i) . (4.10d) 
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Note that ITq is determined with a better precision than IT? (but we shall not need this 
higher precision and give it for completeness). For the components of J- a , we find 
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(4.11a) 
(4.11b) 



The n^'s and JF a 's can now be substituted into the integrals on both sides of ([4.91). This 
is correct because the support of the integrals is the compact support of the source, which 
is, for a slowly- moving source, entirely located within the source's near-zone Di, where the 
post-Newtonian expansion is valid. Straightforward computations permit us to re-express 
the right-hand-sides of ( f4.9| ) into the form of total time-derivatives. We do not detail here 
this computation which is well-known (at 1PN order), but we present in Sect. IV. B a some- 
what general formula which can be used to reach elegantly the result (see Eq. (4.23)). By 
transfering the total time-derivatives to the left-hand-sides of ( |4.9|) , one obtains the looked- 
for conservation laws at 1PN order, namely 
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where the 1PN energy E 1 ™, linear momentum P/ PN , and angular momentum 5*/ PN are 
given by the integrals over the source 
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The 1PN energy E 1FN can also be written as 6^ 
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A similar but more precise computation would yield the 2PN-conserved quantities |27 



B. Secular losses of the lPN-accurate energy and momenta 



As the reactive potentials V^ reac manifestly change sign in a time reversal, they are ex- 
pected to yield dissipative effects in the dynamics of the system, i.e. secular losses of its 
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total energy, angular momentum and linear momentum. The "Newtonian" radiation reac- 
tion force is known to extract energy in the system at the same rate as given by the Einstein 
quadrupole formula, both in the case of weakly self-gravitating systems p3|-|3"9|j and com- 
pact binary systems [[4"0"l-[42|j . Similarly the reaction force extracts angular momentum in the 
system. As concerns linear momentum the Newtonian reaction force is not precise enough, 
and one needs to go to 1PN order. 

In this subsection, we prove that the lPN-accurate reactive potentials l/^ eac lead to 
decreases of the lPN-accurate energy and momenta (computed in ( |4.15|) -( f4.18| )) which are 
in perfect agreement with the corresponding far-zone fluxes, known from the works [B3",4"9",51J 



in the case of the energy and angular momentum, and from the works |55|- p7| , ^9| in the case 
of linear momentum. 

We start again from the equations of motion ( (4.6| )-( p~8|) , which imply, after spatial inte- 
gration, the laws (fOl) that we recopy here: 
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(4.19b) 
(4.19c) 



The left-hand-sides are in the form of total time-derivatives. To 1PN order, we have seen 
that the right-hand-sides can be transformed into total time-derivatives, which combine 
with the left-hand-sides to give the lPN-conserved energy and momenta. Here we shall 
prove that the contributions due to the reactive potentials in the right-hand-sides cannot 
be transformed entirely into total time-derivatives, and that the remaining terms yield pre- 
cisely the corresponding 1PN fluxes. The balance equations then follow (modulo a slight 
assumption and a general argument stated below). 

The right- hand- sides of (|4.19| ) are evaluated by substituting the metric (3.1), involving 
the potentials V™ = V* n + V^ eac . The components of the force density (O) are found to be 
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(4.20a) 
(4.20b) 



where we have been careful at handling correctly the un-controlled 2PN and 3PN approxi- 
mations, which lead to the terms symbolized by the c~ n n J r ^s. The equations ( f4.20p reduce 
to fl4.11|) at the 1PN approximation. They give the components of the force as linear func- 
tional of V m and V™. The remainders are 4PN at least. The same computation using the 
stress-energy density ( |4.7| ) yields the term which is further needed in ( |4.19c| ), 
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The n7?s represent the 2PN and 3PN approximations. Thanks to ( |4.20|) -( |4~2T 
transform the laws (|4.19|) into 



one can now 
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it {- / An s} = / A {-° 8 > v " + ^ d ' y r] + + + (?) • 
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+^ + ^ + o(l) , (4.22c) 
where n X, n Yj, and ra Zj denote some spatial integrals of the 2PN and 3PN terms in ( |4.20| )- 

O). 

Consider first the piece in V' n which is composed of the potential V^ n , given by the 
symmetric integral ( |3.4| ) or by the Taylor expansion ( |3.5| ). To 1PN order V* n contributes to 
the laws ( 4.22|) only in the form of total time-derivatives (see Sect IV. A). Here we present a 



more general proof of this result, valid formally up to any post-Newtonian order. This proof 
shows that the result is due to the very structure of the symmetric potential V™ as given 
by ( p.5| ). The technical formula sustaining the proof is 

a f N—i ) 

o,{x)d?a v {x') + (-) N+1 a v (x')d?a„(x) = | j £ {-fdja,{x)df """V^x') \ , (4.23) 

where x = (x, t) and x' = (x', t) denote two field points (located in the same hypersurface 
t = x°/c = const), and where N is some integer [we recall the notation o~ M = (cr, c^)]. The 
contributions of the symmetric potential in the energy and linear momentum laws (|4.22a|) 
and fl4.22bp are all of the same type, involving the spatial integral of a^{x)d a V^{x) (with no 
summation on \i and a = 0,i). One replaces into this spatial integral the potential V^ n by its 
Taylor expansion when c — > oo as given by ( |3.5| ) . This yields a series of terms involving when 
the index a = the double spatial integral of |x — x / | 2p_1 cr M (x)9 t 2p+1 cr M (x'), and when a = i 
the double integral of <9j|x — ~x!\ 2p ~ l a ll {x)d^ p a l j i {x'). By symmetrizing the integrand under the 
exchange x «-> x' and by using the formula ( |4.23j ) where N = 2p + 1 when a = and N = 2p 
when a = 1, 2, 3, one finds that the integral is indeed a total time-derivative. The same is 
true for the symmetric contributions in the angular momentum law ( 4.22cj ), which yield a 



series of integrals of EijkXjdk\^ — x.'\ 2p ~ 1 a fJ ,(x)d 2p a II (x') and %jt|x — x.'\ 2p ~ 1 aj(x)dt P ak{x'), on 
which one uses the formula ( |4.23j ) where N = 2p. 

Thus the symmetric (inner) potentials V^ n contribute to the right-hand-sides of (f4.22|) 



only in the form of total time-derivatives. This is true even though, as noticed earlier, the 
potentials V^ n contain some reactive ("time-odd") terms, through the contributions of the 
source densities a M . As shown here, these time-odd terms combine with the other time-odd 



terms present in the <t m 's appearing explicitly in ( f4.22|) to form time-derivatives. Such time 



odd terms will not participate ultimately to the balance equations, but they do participate 
to the complete 3.5PN approximation in the equations of motion of the system. This fact 
has been noticed, and these time-odd terms computed for binary systems, by Iyer and Will 
H (see their equations (Of) and (|3.9|)). 
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The numerous time-derivatives resulting from the symmetric potentials are then trans- 
ferred to the left-hand-sides of ( 4.22j ). To 1PN order these time-derivatives permit recon- 
structing the lPN-conserved energy and momenta E 1PN , p. 1PN and S , / PN . We include also 
the time-derivatives of higher order but they will be negligible in the balance equations (see 
below). Therefore, we have proved that the laws ( |4.22| ) can be re- written as 
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(4.24a) 



(4.24b) 



(4.24c) 



The O-symbols 0(c -4 ) denote the terms, coming in particular from the symmetric potentials 
in the right-hand-sides, which are of higher order than 1PN. We add an overbar on these 
remainder terms to distinguish them from other terms introduced below. 

Now recall that the 2PN and 3PN approximations, including in particular the terms 
n X, n Yi, and n Zi in Q4.24|) , are non-radiative (non-dissipative). Indeed they correspond to 
"even" approximations, and depend instantaneously on the parameters of the source. In the 



case of the 2PN approximation, Chandrasekhar and Nutku p7 have proved explicitly that 



±X, 4Y1, and can be transformed into total time-derivatives, leading to the expressions 
of the 2PN-conserved energy and momenta. Here we shall assume that the same property 
holds for the 3PN approximation, namely that the terms qX, G Yi, and §Zi can also be 
transformed into time-derivatives. This assumption is almost certainly correct. The 3PN 
approximation is not expected to yield any secular decrease of quasi-conserved quantities. 
It can be argued, in fact, that the 3PN approximation is the last approximation which is 
purely non-dissipative. Under this (slight) assumption we can now transfer the terms n X, 
n Yi, and n Zi to the left-hand-sides, where they modify the remainder terms 0(c -4 ). Thus, 
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(4.25a) 
(4.25b) 

(4.25c) 



1PN 



+ 



where 0(c 4 ) denotes the modified remainder terms, which satisfy, for instance, E 
0(c" 4 ) = E 2PN + 0(c" 5 ). 

The equations ( |4.25|) clarify the way the losses of energy and momenta are driven by 
the radiation reaction potentials. However, these equations are still to be transformed using 
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the explicit expressions ( |3.6|) - (|378|) . When inserting these expressions into the right-hand- 
sides of ( 4.25|) one is left with numerous terms. All these terms have to be transformed 
and combined together modulo total time-derivatives. Thus, numerous operations by parts 
on the time variable are performed [i.e. Ad t B = d t (AB) — Bd t A], thereby producing many 
time-derivatives which are transferred as before to the left-hand-sides of the equations, where 
they modify the 0(c _4 )'s by some contributions of order c -5 at least (since this is the order 
of the reactive terms). During the transformation of the laws Q4.25a| ) and ( |4.25c| ) for the 
energy and angular momentum, it is crucial to recognize among the terms the expression 



of the lPN-accurate mass quadrupole moment 7y given by ( |3.7| ) (or ( |3.21a| ) with 1 = 2), 
namely 



(4.26) 



And during the transformation of the law (|4.25b[) for linear momentum, the important point 
is to remember that the lPN-accurate mass dipole moment Jj, whose second time- derivative 
is zero as a consequence of the equations of motion [d 2 Ii/dt 2 = 0(c~ 4 )], reads 
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(4.27a) 



This moment is also a particular case, when I = 1, of the general formula (|3.21a|) . An 
alternative expression of the dipole moment is 
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(4.27b) 



which involves the Newtonian conserved mass density (given also by the three first terms in 
(4.15)). Finally, the end results of the computations are some laws involving in the right- 
hand-sides some quadratic products of derivatives of multipole moments, most importantly 
the IPN quadrupole moment ( |4.26| ) (the other moments being Newtonian), namely 
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(4.28a) 
(4.28b) 



+ 



(4.28c) 



The remainders in the left-hand-sides are such that 0(c -4 ) = 0(c~ 4 ) +0(c -5 ). The remain- 
ders in the right-hand-sides are 0(c~ 8 ) in the cases of energy and angular momentum because 
of tail contributions (see Sect. IV. C), but is 0(c~ 9 ) in the case of the linear momentum. 

The last step is to argue that the unknown terms in the left-hand-sides, namely the total 
time-derivatives of the remainders 0(c -4 ), are negligible as compared to the controlled terms 
in the right-hand-sides, despite their larger formal post-Newtonian order (c -4 vs c -5 and 
c -7 ). When computing, for instance, the time evolution of the orbital phase of inspiralling 
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compact binaries P-plj], one uses in the left-hand-side of the balance equation the energy 
valid at the same post-Newtonian order as the energy flux in the right-hand-side. Because 
the difference between the orders of magnitude of the two sides of the equations is c~ 5 , we 
need to show that the time-derivative increases the formal post-Newtonian order by a factor 
c -5 . In ( pT27D this means dO(c~ A )/dt = 0(c' 9 ) [actually, 0(c~ 8 ) would be sufficient in (4.28), 
but 0(c~ 9 ) will be necessary in Sect. IV. C]. In the case of inspiralling compact binaries, such 
an equation is clearly true, because the terms 0(c -4 ) depend only on the orbital separation 
between the two bodies (the orbit being circular), and thus depend only on the energy 
which is conserved at 2PN order (for non-circular orbits one would have also a dependence 
on the angular momentum). Thus the time-derivative adds, by the law of composition of 
derivatives, an extra factor c -5 coming from the time-derivative of the energy itself. More 
generally, this would be true for any system whose 2PN dynamics can be parametrized 
by the 2PN-conserved energy and angular momentum. This argument could perhaps be 
extended to systems whose 2PN dynamics is integrable, in the sense that the solutions 
are parametrized by some finite set of integrals of motion, including the integral of energy. 
Another argument, which is often presented (see e.g. [0), is that the terms dO (c~ 4 ) /dt are 
negligible when taken in average for quasi-periodic systems, for instance a binary system 
moving on a quasi-Keplerian orbit. The time average of a total time-derivative is clearly 
numerical small for such systems, but it seems difficult to quantify precisely the gain in order 
of magnitude which is achieved in this way, for general systems. The most general argument, 
valid for any system, is that the terms dO(c~ A ) / dt are numerically small when one looks at 
the evolution of the system over long time scales, for instance At ^> 0(c~ 4 )(dE 1PN /dt)" 1 
(see Thorne [jBlfl , p. 46). 

Adopting here dO(c~ A )/dt = 0(c~ 9 ) and the latter general argument, we can neglect the 
terms 0(c -4 ) and arrive to the 1PN energy-momenta balance equations 
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relating the lPN-conserved energy and momenta, given by the explicit integrals over the 
source fl3.17| )-( p.l9D , to some combinations of derivatives of multipole moments, also given 



by explicit integrals over the source (see (|3.7|)-(|3~8|)). Note that at this order both sides of 
the equations are in the form of compact-support integrals. The right-hand-sides of (|4.29Q - 
( |1.31| ) agree exactly with (minus) the fluxes of energy and momenta as computed in the 



wave zone of the system. See for instance the equations (4.16'), (4.20') and (4.23') in [49 
when truncated to 1PN order [and recalling that the quadrupole moment which enters the 
1PN fluxes is precisely the one given by (4.26)]. Thus, we can conclude on the validity of 
the balance equations at 1PN order, for weakly self-gravitating systems. 



These equations could also be recovered, in principle, from the relations (|2.14fj -( p.l5 ) 



(which were obtained in paper I). Indeed (2.14)-(2.15) involve, besides some instantaneous 
contributions such as T^(t), some non-local (or hereditary) contributions contained in the 
functions m(t), rrii(t) and Sj(i). These contributions modify the constant monopole and 
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dipole moments M, Mi and Si by some expressions which correspond exactly to the emitted 
fluxes. The balance equations could be recovered (with, though, less precision than obtained 
in this paper) by using the constancy of the monopole and dipoles M, Mj and Si in the 
equations (2.14) written for I = and I = 1, and by using the matching equations obtained 
in ( |3.27| ), also written for I = and 1 = 1. Related to this, notice the term involving a single 
time-antiderivative in the function mi(t) of (|2.15b| ), and which is associated with a secular 
displacement of the center of mass position. 



C. Tail effects at 1.5PN order 



To 1.5PN order in the radiation reaction force appears an hereditary integral (i.e. an 
integral extending on the whole past history of the source), which is associated physically 
with the effects of gravitational- wave tails. More precisely, it is shown in [53~|, using the 
same combination of approximation methods as used in paper I and this paper, that the 
dominant hereditary contribution in the inner post-Newtonian metric (valid all over Di) 
arises at the 4PN order. At this order, the dynamics of a self-gravitating system is thus 
intrinsically dependent on the full past evolution of the system. 

In a particular gauge (defined in p3"|), the 4PN-hereditary contribution in g™ u is entirely 
located in the 00 component of the metric, and reads 
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(4.32) 



The hereditary contributions in the other components of the metric (Oi and ij) arise at higher 



order. Note that the hereditary (tail) integral in ( 4.32 ) involves a logarithmic kernel. A 
priori, one should include in the logarithm a constant time scale P in order to adimensionalize 
the integration variable A, say ln(A/2P). However, InP would actually be in factor of an 
instantaneous term [depending only on the current instant t through the sixth time-derivative 



I-j(t)], so ( |4.32| ) is in fact independent of the choice of time scale. In (|4.32|) we have chosen 
for simplicity P = 1 sec. The presence of the tail integral ( |4.32j ) in the metric implies a 
modification of the radiation reaction force at the relative 1.5PN order 53fl. The other 4PN 



terms are not controlled at this stage, but are instantaneous and thus do not yield any 
radiation reaction effects (indeed the 4PN approximation is "even" in the post-Newtonian 
sense). It was further shown that the 1.5PN tail integral in the radiation reaction is 
such that there is exact energy balance with a corresponding integral present in the far-zone 
flux. Here we recover this fact and add it up to the results obtained previously. 

As the gauge transformation yielding ( f4.32| ) in deals only with 4PN terms, it can be 
applied to the inner metric g™ u given by ( |3.1|) without modifying any of the known terms at 
the 1PN non-radiative and reactive approximations. It is clear from ( [4.32 ) and the reactive 
potentials (|3.6| ) that after gauge transformation, the inner metric takes the same form as 
(ft.l|), except that the reactive potentials are now more accurate, and given by 
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Still there remain in the metric some un-controlled (even) 4PN terms, but these are made 
of instantaneous spatial integrals over the source variables, exactly like the un-controlled 
2PN and 3PN terms. [The expressions ( |4.33| ) can be recovered also from Sect. III.D of 
paper I and a matching similar to the one performed in this paper.] With ( |4.33| ) in hands, 
one readily extends the balance equations to 1.5PN order. First one obtains (4.24), but 
where the reactive potentials are given more accurately by ( 4.33 ), and where there are some 
instantaneous 4PN terms sX, sYi and 8%i in the right-hand-sides. Extending the (slight) 
assumption made before concerning the similar 3PN terms, we can transform 8 X, 8 Y^ and 
8 Zj into time-derivatives and transfer them to the left-hand-sides. This yields (4.25), except 
that the remainders in the right-hand-sides are 0(c~ 9 ) instead of 0(c~ 8 ). Using Q4.33D , 
we then obtain (working modulo total time-derivatives) the laws (4.28) augmented by the 
tail contributions arising at order c -8 in the right-hand-sides. The remainders in the left- 
hand-sides are of the order dO(c~ A ) / dt = 0(c -9 ) (arguing as previously), and therefore are 
negligible as compared to the tail contributions at c -8 . In the case of energy the 1.5PN 
balance equation is obtained as 



dE 



1PN 



dt 



G r(3)r(3) 
5c 5 H ij 

4G 2 M 



1 




r(4) , t(3) t(3) 

1 ijk 1 ijk ^ J ij J ij 



(4.34) 



Because there are no terms of order c -3 in the internal energy of the system (see Sect. IV. A), 
the energy E lPN appearing in the left-hand-side is in fact valid at the 1.5PN order. Finally, 
to the required order, one can re- write ( (4.34 ) equivalently in a form where the flux of energy 
is manifestly positive-definite, 



dE 



1PN 



dt 



2GM 



(T<A\2 + ™(jQh 

lg9 V ijk) T 45 ^ij , 




d\ In 



+ 




A) 



(4.35) 



Under the latter form one recognizes in the right-hand-side the known energy flux at 1.5PN 
order. Indeed the effective quadrupole moment which appears in the parenthesis agrees with 
the tail-modified radiative quadrupole moment parametrizing the field in the far zone (see 
Eq. (3.10) in fJ5|). [The term associated with the (gauge-dependent) constant 11/12 in the 
radiative quadrupole moment [45| yields a total time-derivative in the energy flux (as would 
yield any time scale P in the logarithm), and can be neglected in (|4.35|) .1 The 1.5PN balance 
equation for angular momentum is proved similarly (it involves as required the same tail- 
modified radiative quadrupole moment). The balance equation for linear momentum does 
not include any tail contribution at 1.5PN order, and simply remains in the form (00). 
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